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Abstract,  The  steady-state  propagation  of  a  semi -infinite,  anti -plane 


shear  crack  is  reconsidered  for  a  general,  infinite,  homogeneous  and 
isotropic  linearly  viscoelastic  body.  As  with  an  earlier  study,  the 
inertial  term  in  the  equation  of  motion  is  retained  and  the  shear  modulus 
is  only  assumed  to  be  positive,  continuous,  decreasing  and  convex.  A 
Barenblatt  type  failure  zone  is  introduced  in  order  to  cancel  the  singular 
stress,  and  a  numerically  convenient  expansion  for  the  dynamic  Energy 
Release  Rate  (ERR)  is  derived  for  a  rather  general  class  of  crack  face 
loadings.  The  ERR  is  shown  to  have  a  complicated  dependence  on  crack 
speed  and  material  properties  with  significant  qualitative  differences 
between  viscoelastic  and  elastic  material.  The  results  are  illustrated 
with  numerical  calculations  for  both  power-law  material  and  a  standard 
linear  solid. 


§1  Introduction 


A  central  issue  in  fracture  mechanics  is  the  development  of  fracture 
criteria.  A  great  many  experimental  and  analytical  studies  have  addressed 
this  topic  in  the  nearly  sixty  years  since  Griffith's  pioneering  work. 

One  fact  that  has  emerged  from  this  effort  is  that  the  choice  of  a 
fracture  criterion  is  very  much  dependent  upon  the  particular  scenario 
considered.  For  example,  the  notion  of  a  critical  Stress  Intensity  Factor 
(SIF)  has  provided  a  highly  successful  criterion  for  quasi-static  crack 
propagation  in  linearly  elastic  material.  Important  factors  for  the 
success  of  the  SIF  in  this  setting  are  that  it  is  often  easily  computed 
and  that  the  Energy  Release  Rate  (ERR)  can  be  determined  in  a  simple 
manner  from  it.  However,  such  is  not  the  case  for  dynamically  propagating 
cracks  in  viscoelastic  material  for  which  the  entire  deformation 
history  during  crack  advance,  not  just  the  SIF,  is  required  in  order  to 
calculate  the  ERR. 

The  notion  of  ERR  employed  here  provides  a  phenomenologically 
meaningful  and  mathematicaly  convenient  fracture  criterion  for  dynamic 
viscoelastic  crack  propagation.  For  dynamic  fracture  in  linearly  elastic 
material,  the  ERR,  G,  is  usually  defined  as  the  difference  between  the 
total  power  input  to  the  body  by  external  forces,  P,  and  the  rate  of 

change  of  the  total  internal  strain  energy,  E,  and  kinetic  energy,  K.  It 
is  then  easily  shown  (see  [1]  for  example)  that  G  can  be  calculated  from  a 
knowledge  of  the  instantaneous,  singular  asymptotic  stress  field  at  each 
crack  tip. 

For  viscoelastic  material  the  situation  is  more  complicated.  In  this 
case  the  difference  P-E-K  contains  two  terms:  a  history  dependent  body 


'sJk 


term,  W,  representing  the  energy  dissipated  through  viscous  effects  and 
the  quantity  G  which  now  depends  upon  the  entire  history  of  the  singular 
asymptotic  stress  and  strain  fields  at  each  crack  tip  during  the  time  the 
tip  is  advancing.  Physically  it  is  reasonable  to  assume  that  a  fracture 
criterion  should  depend  only  upon  these  near  tip  fields.  However,  the 
quantity  G,  which  can  be  interpreted  as  the  power  input  to  the  crack  tip, 
is  extremely  cumberson  to  calculate  from  the  singular  fields.  Instead,  in 
this  paper  a  Barenblatt  type  failure  or  process  zone  is  introduced  behind 
the  crack  tip  in  order  to  cancel  the  singular  asymptotic  fields  in  front 
of  the  crack.  What  results  is  that  G  is  given  by  the  formula  (3)  below 
and  hence  that  the  power  input  to  the  crack  tip  is  given  by  the  rate  of 
work  done  by  the  tractions  in  the  failure  lone. 

The  purpose  of  this  paper  is  to  demonstrate  that  for  dynamic  crack 
propagation  in  linear  viscoelastic  material,  the  ERR  can  be  conveniently 
calculated  even  though  it  is  not  just  a  simple  function  of  the  SIF. 
Moreover,  while  the  SIF  is  a  monotomical ly  decreasing  function  of  crack 
speed,  the  ERR  can  exhibit  much  more  complicated  behavior  depending  upon 
combined  viscoelastic  and  dynamic  effects.  In  particular,  the  numerical 
calculations  presented  here  suggest  that  for  some  ranges  of  crack  speed, 
stable  crack  growth  occurs  while  for  other  values,  unstable  propagation 
can  be  expected.  This  point  is  considered  further  when  the  numerical 
examples  are  discussed. 


The  specific  boundary  value  problem  considered  here  is  that 
corresponding  to  the  steady  propagation  (to  the  right)  with  speed  V  of  a 
semi -infinite,  anti-plane  shear  crack  in  a  general,  homogeneous  and 
isotropic,  linearly  viscoelatic  body.  The  shear  modulus,  ii(t),  is  assumed 
only  to  be  a  positive,  non -increasing  and  convex  function  of  time,  t.  The 
governing  equation  of  motion  for  the  out  of  plane  displacement,  u^,  is 

u*dAU3  =  p  U3^^^ 

2  2  2  2 

where  A  is  the  two-dimensional  Laplacian,  A  =  (9  /9Xj^)  +  (9  /9X2),  and 

u*  de  denotes  the  Riemann-Stieltjes  convolution 

t 

u*  de  =  /  u(t-T )  de(T )  . 

Upon  adoption  of  the  Galilean  variables  x  =  Xj  -  Vt,  y  «  Xg,  the  boundary 
conditions  may  be  written 

a23(x,0)  =  ly  (w*  dUa)  =  f(x),  x<0 
U3(x,0)  =  0  ,  X  >  0 

2  2 

aij(x,y)  -*-0,  X  +  y  -►  » 

where  a  . .  are  the  stress  components  and  f (x)  is  a  system  of  tractions 
^  J 

moving  with  the  crack. 

The  starting  point  of  the  present  investigation  is  the  solution 

derived  in  [3]  for  the  above  boundary  value  problem.  It  was  shown  in  [3] 

★ 

that  two  cases  arise  naturally  in  constructing  the  solution:  0  £  V  <  C 

★  *  1 /7  1/2 

and  C  <  V  <  C  where  C  =  (u(«)/p)  and  C  =  (u(0)/p)  are  the  elastic 

shear  wave  speeds  corresponding  to  the  equilibrium  and  glassy  values  of 


the  shear  modules  u(t).  For  0  £  V  <  C  ,  the  stress  field  is  that  for 
static  elastic  fracture  and  is  therefore  independent  of  crack  speed  and 

•k 

material  properties.  Whereas,  for  C  <  V  <  C,  the  stress  field  is  both 
speed  and  material  dependent.  Specifically,  the  SIF,  K,  is  given  by 

0  zL 

I  023'^^*°^  |x|  ^  dx  ,  0  <  V  <  C*  ( 

1  ^  ^  * 

-  i  /  a23'(x,0)  |xr  e  °  dx  C  <  V  <  C 

«O0 

where  g"(x)  (g^(x))  denotes  the  restriction  of  g(x)  to  x  <  0  (x  <  0)  and 
is  the  unique  positive  constant  such  that 

Vq-  /  P  (t)e  °  dt  =  {\l/cr  .  (: 

°  0 

k 

In  (2),  y  (t)  denotes  the  nondimensional  modulus  given  by 

u  (t)  =*  ij{t)/ii(0) . 

In  order  to  calculate  the  ERR,  the  object  of  principal  interest  in 
the  present  study,  it  is  necessary,  as  discussed  earlier,  to  modify  the 
above  boundary  value  problem  by  the  introduction  of  a  Barenblatt  type 
failure  zone.  Specifically,  it  is  now  assumed  that  two  loads  are  acting 

on  the  crack  faces:  the  applied  (external)  tractions  a23(x,0)  discussed 
above,  but  now  denoted  ag(x),  and  cohesive  (failure)  stresses 
a^(x)  acting  in  a  failure  zone  of  length  a^  immediately  behind  the  crack 
tip.  The  only  assumptions  about  o^(x)  are  that  a^  is  small  relative 
relative  to  some  length  scale  a^  associated  with  o~(x)  and  that 


Kg  +  Kf  =  0  where  and  are  the  SIF's  corresponding  to 

Og  and  Op  respectively.  Hence  the  effect  of  the  failure  zone  is  to 

cancel  the  singular  stresses  ahead  of  the  crack  tip  and  thereby  produce 
a  cusp  shaped  crack  profile  behind  the  tip. 

The  ERR,  G,  for  steady-state  crack  propagation  is  now  defined  by 

0 

G  =  /  a~(x)  u"  j(x,0)dx  (3) 

“^f 

where  U2(x,0)  is  the  crack  face  displacement  corresponding  to  the  combined 

loading  Og  +  •  Thus  G  has  the  interpretation  of  power  input  to  the 

crack  tip  from  the  total  load,  i  .e.  the  power  available  to  the  crack 
tip  for  propagating  the  crack.  Unlike  for  elastic  material,  G  for 
viscoelastic  material  is  not  merely  a  simple  function  of  K.  Rather,  as  is 

evident  from  [3],  u^  j(x,0)  has  a  complicated  dependence  upon  the  loading 

ag(x,0),  making  impracticable  the  direct  numerical  evaluation  of  (3).  In 

the  next  section,  a  computationally  convenient  expression  for  (3)  is 
derived  for  a  special,  but  still  fairly  general,  class  of  loadings 

o”  and  0^. 

§2.  The  Calculation  of  G 

For  simplicity  of  argument  and  clarity  of  result,  the  ERR,  G,  given 
in  (3)  will  be  calculated  first  for  a  simple  special  case.  Specifically, 

the  external  load,  a”(x)  and  failure  zone  stresses,  a^(x)  will  be  assumed 


to  have  the  forms 


6 


x/a 

Og(x)  =  Lg  e 

-  «  <  X  <  0  (4) 

x/a- 

af(x)  =  -  e 

where  a^/a^  «  1.  (c.f.  [2])  For  a^r/ag  small  enough,  the  fact  that  (x) 

does  not  have  support  in  some  small,  compact  interval  behind  the  crack  tip 
will  have  a  negligible  effect  on  the  results.  The  assumptions  (4)  then 
clearly  incorporate  the  salient  features  of  the  Barenblatt  model,  namely, 
a  set  of  cohesive  stresses  and  associated  length  scale  a^  and  a  length 

scale  a  associated  with  the  applied  load  a"  such  that  al  cancels  the 
0  0  1 

singular  stresses  produced  bya~  and  a^/a^  «  1.  If  should  be  noted  that 

in  this  case,  (3)  is  replaced  by 

0 

G  =  J  a^(x)  U3^j(x,0)dx.  (5) 

In  order  to  present  the  results  in  a  suitable  nondimensional  form,  it 
is  useful  to  introduce  certain  parameters.  First  a  characteristic  time 
scale,  T,  is  defined  and  the  shear  modulus  given  the  form 

u(t)  =  m(t/T  ) 

where  m(s)  is  a  nondimensional  function  of  s  with  lim  m(s)  =  1  and 

=  lim  M(t).  Also  useful  are  the  nondimensional  parameters  a,  $,  y,  e 

t-H» 

and  A  defined  by 

a  =  C  x/a^  B  =  q^a^  y  =  V/C  e  =  a^/a^  A  =  L^/L^.  (6) 

It  should  be  nr>ted  that  a  and  6  assume  any  positive  value  while  y  must  be 

★  ★ 
such  that  0  <  Y  <  C/C  .  Also,  B  =  0  whenever  0  <  V  £  C  . 
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The  desired  expression  for  G  is  most  easily  expressed  in  terms  of  the 


Carson  transform,  m(s),  of  m(t)  defined  by 


m(s)  =  m(0)  +  /  e“^^  dm(t). 

0 

It  is  useful  to  record  for  future  reference  the  easily  derived  formula 


p(r)  =  li^m(rT). 

It  will  now  be  shown  that  for  loads  of  the  form  (4),  G  is  given  by 


r  =  1  r 

^  ^  m(aY/e)  L 1-Y^/ml 


.2  T  2 


m(aY/e)  Ll-T  /m(aY/e)  J 

where  K  =  =  -K^  is  the  dynamic  SIF.  From  (1)  and  (4)  it  follows  easily 


K  .  =  '-JIS.  . 

/l+ge  /T+B 

A 

The  derivation  of  (7)  from  (5)  utilizes  the  Fourier  transform  f 
of  a  function  f.  Specifically, 

f(p)  =  /  e’^^dx 


with  inverse,  f(x),  given  by 


f(x)  =  -^  /”  ^'(P)  e'^’^'^dp. 

.00 

Applying  the  Parseval  formula  for  the  Fourier  transform,  it  follows  from 
(5)  that 

00  A 

G  =  /  o"(p)u3^j(p)dp.  {< 
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In  consideration  of  (4),  a  straight  forward  calculation  shows  that 


af(p)  = 


(1+ia^p) 


-afLf 


®f^P^  "  2ir(l-ia^p)  • 

The  integral  in  (9)  may  be  readily  evaluated  using  residues  since,  from 
(10),  it  is  clear  that  a^(p)  has  a  meromorphic  extension  to  the  lower 
complex  half-plane  with  a  simple  pole  at  -i/a^.  Moreover,  u^  2^(x,0) 

vanishes  for  x  >  0  from  which  it  follows  that  03^1  (p)  =  U3j(p)  has  an 
analytic  extension,  F~(z),  to  the  lower  half-plane  with 


lim  F"(p-iq)  =  0. 


Consequently,  for  6  one  has 


G  =  L^F“(-i/a^).  (1 

It  remains  to  evaluate  F'(-i/a^).  To  this  end  use  will  be  made  of 
the  following  formulas  derived  in  [3]: 

U3^l(p)  =  o{p)/G(p),  (1 

G(p)  =  -  i  sgn(p)  ■i7(iVp)Yj(iVp) ,  (] 


Yi(ipV)  =  (l-Y^/u(iVp))^^^ 


Here  one  has 


a  (p)  =  a'(p)  +  o^(p) 


0  '(x)  =  oe(x)  +  of(x) 


where  o^(x)  (a"(x))  denotes  the  restriction  ofo(x)  to  the  half-line 

X  >  0  (x  <  0).  It  should  be  noted  that  o^(x)  is  the  nonsingular  stress 

field  ahead  of  the  crack  that  results  from  superposing  a"  and 

Moreover,  it  is  shown  in  [3]  that 

a''^(p)  =  lim  F''’(z) 

l„(z).0* 

where 


F'*’(z)  =  (q  -iz)^^^  /  o‘(T)(q  -it) 

°  2iri  -« 


-1/2  dr 


(t-z) 


(15) 


with  q^  given  in  (2).  To  evaluate  (15),  consider  first 


a!(p)  =  lim  F^(z)  , 


I„(z)-0+ 
m 


in  which 


~  /  ol  (T){q--iT)"^^^  .  (16) 


2iri  -• 


(t-z) 


From  the  analog  of  (10)  forOg(p),  one  sees  that  o”(p)  has  a  simple  pole 


-1/2 

at  i/a  .  Since  the  branch  of  (q  -it)  '  must  be  chosen  to  be  analytic  in 
e  0 

the  upper  half-plane  (see  [3]),  the  integral  in  (16)  can  be  calculated 
using  residues.  Thus  for  I^(z)  >  0, 

Fg(z)  =  -(qo-iz)^^^[-iLg(qQ+l/ag)'^''^/(i/ag-z)  +  a‘(z)(qQ-iz)‘^/^]. 

Now  letting  I^(z)  -►  0  and  adding  o“(p)  too  (p)  there  results 


Oe(p)  =  F'^(p)  +  o"(p) 

=  i  Lg(qQ-ip)^^^(qo+l/ag)'^^^/(i/ag-p) 

=  agLg(6/(l+0))^^^(l-ip/qQ)^^^/(l+ipag). 


(17) 


Similarly,  foro^(p)  one  can  show  that 


Oflp)  =  -a^L^(Be/(l+ee))^''^(l-i  p/q^)^^^/(l+ipa^) . 


Combining  (4),  (8),  (17)  and  (18)  one  concludes  that 


a(p)  =  -a  L  (e/(i+e))^/^ 


(l-e)(l-ip/q^)^/^ 


(l+iagp)(l+ia^p) 


Substitution  of  (19)  into  (12)  gives 


^  ^  agLg(3/(l+0))^^^  (l-e)*j(p)*2(p) 


*j(p)  =  (l+iagp)‘^(l+ia^p)‘^ 


♦2(P)  =  (l-ip/qQ)^'VG(p). 

't^(p)  is  obviously  an  analytic  function  in  the  lower  half-plane  and  in  [3] 


it  was  shown  that  the  branches  of  (l-ip/q^)^'^  and  Yj(iVp)  must  be  chosen 
so  that  analytic  there.  Thus  one  may  substitute  (20)  into 


(11)  making  use  of  (13)  and  (14)  to  conclude  that 


1/2/, 


l-l/(q„aJ  1  1/2 


2(l+ag/a^)u(V/a^)  L  1-Y  /m(V/a^) 


Combining  (6),  (8),  and  (21),  there  finally  results  the  desired  formula 


The  derivation  given  above  is  easily  modified  to  produce  an  analog  to 


formula  (7)  for  more  general  loads  of  the  form 


V.  <.  •*, 


=  L  /  e 
e  e  Q 


d  hg(t) 


«  tx/a- 

<i'Ax)  =  -L,  /  e  ^  d  h.(t) 

'  ^0  ' 

where  hg(t)  and  h^(t)  are  arbitrary  signed  (not  necessarily  positive) 

measures  restricted  only  to  the  extent  that  the  required  integrals 
converge.  For  example,  the  special  cases  (4)  correspond  to  d  hg{t)  = 

d  h^(t)  =  6(t-l),  the  Dirac  measure  concentrated  at  t  =  1.  As  further 

examples,  d  hg(t)  =  sin(t)dt  and  d  h^(t)  =  cos(t)dt  produce 

1 

= - 2  • 

IMx/a^)"^  *  l+(x/a^r 

respectively. 

Substitution  of  (22)  into  (1)  followed  by  an  interchange  of  order  of 
integration  results  in 


=  Lg/a^  /  {t+e)‘  d  h^(t) 


"  ‘■f’^  I  d  h^(t) 


=  -  K. 


where,  as  before,  6  =  a^q^  .  Corresponding  to  (10)  there  is 


!(P)  =  ^pLp  /  (t+ipa^)  d  hg(t) 


o'jp)  -  4~/ 


•/  (t-ia  o)’  d  h^(t)  . 


^  .•  j'  j-  .-  J  -r  ^  i-  ^  <  j  X O  •J.O  •/  •.'••V.'.V.Vi.'.' 
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A 

Lines  (15)  and  (16)  are  still  valid  but  with  given  now  by  (23).  In 
particular,  one  has  for  Ip,(z)  >  0  and  after  an  interchange  of  integration 
that 

Fe(2>  '  ST  .  (24) 

The  integrand  in  the  inner  integral  in  (24)  is  analytic  for  I  (t)  >  0 

m 

except  for  simple  poles  at  t  =  i  t/a^  and  t  =  z.  Calculating  the  inner 
integral  by  residues,  there  then  results 

If  one  now  lets  I^(z)  +  0  in  this  last  result  and  makes  use  of  (23a)  and 
(24)  it  follows  that 

%(p)  =  "e  +  %  '  4eL^(l-fP/q„)*'^  /  (B/('*6))‘''2  (,4pa^)  •  (25) 
For  o ^(p)  one  has 


^  1  /9  **  1  /9  ^h^(T  ) 

Of(p)  =  -a^L^d-ip/q^)^'^  /  . 

the  analog  of  (18),  which  when  combined  with  (25),  gives 

A  A  A 

O  (p)  =  Og(p)  +  Of(p) 

=  (1-i  P/qo)^^^H(p)  (26) 


with  H(p)  = 


-  a  L  (S/(t+6))'^^  «  a,L,(e6/(t+E8))^(^ 

i  CxFia^p) - ^".(F)-/  ■■  (,.(a,p) - PP,(').(22) 


It  is  easily  seen  that  H(p)  is  analytic  for  I^(p)  <  0  and  H(p)  +  0  as 

I^(p)  After  substitution  of  (11),  (12),  (23b),  (26),  and  (27)  into 

(9)  and  an  interchange  of  integration  one  obtains 

«0  00 

G  =  a^  /  dh^(t)  ^  /  'r(t,p)dp  (28) 

0 

where 


H(p)(l-ip/qn) 
=  (t-ipa7)6(p; 


As  before,  the  branches  of  Yj(ipV)  and  /l-ip/q^  can  be  chosen  so  that 

i  sgn(p)  /rT~p7^/Y^(ipV)  is  analytic  for  I|^(p)  <  0.  Thus,  t(t,p)  is 
analytic  for  Ij„(p)  <  0,  except  for  a  simple  pole  at  p  *  -  it/a^. 

Evaluting  the  inner  integral  in  (28)  by  residues,  yields  finally 

U  -  H(-it/a.)  r  1. 1  1/2 

G=^/  dh^(t)  - 1-  (29) 

°  .  l-Y^/m(aYt/e)_ 

with  H(-it/a^)  defined  in  (27). 

The  two  integrations  required  to  evaluate  (29)  make  it  much  more 
cumbersome  to  calculate  numerically  than  (7),  though  still  much  easier 
than  calculating  G  directly  from  either  (3)  or  (5).  However,  since 
e  «  1,  it  is  not  unreasonable  to  take  dh^(t)  =  6(t-l),  i.e.  to  take  the 

simple  form  (4b)  foro^(x),  since  the  details  of  the  failure  zone  stress 
are  not  significant.  Formula  (29)  then  simplifies  greatly  to 


_  '-f  r  l-l/(ge] 


m(aY/e) 


1-Y  /m(aY/e) 


"3fLf  1/0  1/0  “ 

H(-i/a-)  =  — ^  (e6/(l+e3))^/^  +  a  L  f  - ? - ^ 

^  ^  ®  ®  0  (T+l/e)(T+3)^' 

subject  to  the  auxiliary  constraint 


I  ( _ L. 

l+0e 


a/2  _  ,  /  a/2  f 


In  the  next  section,  the  qualitative  behavior  of  G  is  investigated  by 
considering  the  special  cases  of  a  power>law  material  and  a  standard 
linear  solid. 


§3.  Numerical  Examples 

The  formula  (7)  will  now  be  applied  to  the  special  cases  of  a 
standard  linear  solid  and  power-law  material.  First  considered  is  the 
standard  linear  solid,  which  is  modelled  by  a  constant  Poisson's  ratio  and 
a  shear  modulus,  M(t),  of  the  form 

u(t)  =  u^(l+ne"^^^) 

=  w  m(t/T). 


It  follows  that  1  +  n  =  (C/C*)^  and  m(s)  is  given  by 


m(s)  = 


l+s(l+^ 


From  (2)  and  (30)  one  easily  shows  that 

q,*— 

It  should  be  noted  that  the  restriction  0  <  V  <  C  corresponds  to 


0  <  <  1  +  n ,  and  moreover,  that  qo  “  0  for  0  £  T  <.  !•  From  (7),  (8), 


(13),  (30),  and  (31),  one  readily  obtains  the  formula 


G  = 


g(o,Y,n  ,e) 


where 


g(a,Y,n,e)  =  (|^) 


JI+oy) 


(1+oy(1-Hi  ))^''^(l-Y^+aY(HTi  -y'")) 
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and 


(32 


g(a.Y,n,e)  =  (|^) 


(e+aY)(l+ee) 

2  a/2 


1/2 


(l+6)(oY(n+l-Y  ))  (e+aY(Hn)) 


1/2 


1  <  Y  <  (l-*n ) 


1/2 


Figure  1  displays  a  normalized  SIF,  k  =  K/(Lg/ag),  which  from  (8)  is 


.-1 


seen  to  be  just  k  =  (1+0)"  .  From  (31),  3  is  seen  to  be  given  by 


3  = 


ir-1) 


aY((l-^  )-Y  ) 


In  Figure  1,  k  is  plotted  atainst  y'  =  y//1-hi  for  a  =  .1,1  and  10, 
Clearly,  k  must  vanish  as  y  approaches  1. 

In  Figure  2,  g(a,Y,n  ,g)  is  plotted  against  log(Y)  for  0  <  y  1  and 

against  log((/nT-l)/(/fr+r-Y))  for  1  <  y  <  *^+1  ,  for  n  =  10,  e  =  .01  and 
a  =  .1,  1,  10,100.  Thus  the  failure  zone  length  is  assumed  to  be 
constant.  For  many  materials  (such  as  rubber)  a  more  realistic 
approximation  is  furnished  by  assuming  a  constant  failure  zone  stress 
level,  L^.  This  is  tantamount  to  holding  X  constant  and  allowing  e  to 

vary.  From  (8)  one  easily  calculates  e  as  a  function  of  X  to  be 


e  =  [A^(3+l)-3]’^ 


Consequently,  one  may  regard  g  as  a  function  of  o,  y,  n  and  X. 

Figure  3  is  the  analog  of  Figure  2  for  X=10,  n=10  and  a=.l,  1,  10. 

It  should  be  noted  that  g  vanishes  as  y  approaches  /HT  when  X  is  held 
constant.  However,  with  e  constant,  g  tends  to  a  non-zero  finite  limit  as 
Y  tends  to  /n+l ,  i.e.  as  V  approaches  the  glassy  shear  wave  speed. 

Indeed,  from  (32b)  it  is  easily  seen  that 


where 


g(a,Y,n,e) 


f  (T^)Cl-A(aY)"r^/^[l+(aY)"r^/2  q  <  y  <  1 


^  ^  l-(£0)" 

Figure  4  is  a  plot  of  the  non-dimensional  SIF  k  =  (1+3)"^  against  y  for 
o=.l,  1,  10  and  n=.3.  Figure  5  shows  g  plotted  against  log(Y)  for 
a=.l,  1,  10,  100,  n=.3,  ea.Ol.  Figure  6  has  a=l,  e=.01,  n=.l,  .3,  .5, 

.7.  The  case  X  constant  is  not  exhibited  here  since  it  results  in  little 
change  from  the  constant  e  calculations. 

Several  comments  on  the  numerical  results  should  be  made.  It  can  be 
observed  for  a  standard  linear  solid  in  Figures  2  and  3  and  may  be  shown 
analytically  for  general  material  that  the  slope  of  the  curve  g  versus  y 
is  discontinuous  for  y  =  1,  i .e.  as  V  passes  through  C*.  A  more  striking 
observation  is  the  loss  of  monotonicity  of  g  as  a  function  of  y  for 
certain  ranges  of  the  parameters.  For  example,  in  Figures  2  and  3,  it  is 
seen  that  for  a  =  .1,  g  is  monotone  decreasing  in  y.  whereas  for  a  =  1,  10 
g  has  a  relative  maximum  on  1  <  y  <  /ThT.  As  seen  from  Figure  6,  for 
power-law  material,  varying  the  exponent  n  also  causes  a  transition  from 
monotonicity  to  having  a  single  relative  maximum.  The  lack  of 
monotonicity  suggests  that  certain  crack  speeds  are  unstable.  In 
particular,  since  G  has  the  interpretation  of  power  input  to  the  crack 
tip,  the  Y-intervals  on  which  g  is  increasing  are  those  on  which  an 
increase  in  crack  speed  produces  an  increase  in  the  power  available  to 
propagate  the  crack.  Evidently  this  lack  of  monotonicity  is  due  to  the 
combined  inertial  and  dynamics  effects  considered  here.  It  is  also  worth 
noting  that  these  results  illustrate  that  in  contrast  to  elastic  material. 
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for  viscoelastic  material,  there  is  no  simple  relationship  between  G  and 
K.  Indeed,  though  K  is  always  a  monotone  decreasing  function  of  V  which 
vanishes  at  the  glassy  shear  wave  speed,  G  need  not  be  monotone  and  need 
not  vanish  at  V  =  c. 
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